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1. Introduction 

The Svetitsky-Yaffe conjecture [jl|] states that the Yang-Mills finite temperature transition in dimen- 
sion d + 1 is described by an effective spin model in d dimensions with short range interactions. 
This relationship is analyzed for a SU (3) YM theory with adjoint Polyakov loop potential and the 
SU (4) YM theory in terms of inverse Monte-Carlo (IMC) methods. Ways to perform IMC are 
provided by demon methods ^ ||]. Specific thermalization effects of these must be discussed to 
obtain reliable results. Taking these effects into account we compare the corresponding results to a 
Schwinger-Dyson approach to IMC [§]. 

2. Effective models for Yang-Mills theories 

We start with the well-known lattice Wilson action 

^^(l-^Retrf/o), P = -fe (2.1) 

and perform a strong coupling expansion (for small j3). Since the resulting 'operators' (Polyakov 
loop monomials) are dimensionless there is no natural ordering scheme. We therefore use a trun- 
cation scheme based on the ordering by powers of j8 which are closely related to the dimension of 
the corresponding group representations and ordering by the distance across which the Polyakov 
loops are coupled. In compact form the strong coupling expansion is given by 

Seff = I £ L4i(i3)ri%A- = E^ (2-2) 

r St v ..M r h...i r (-1 i 

with the basic building blocks 

Ssi,i = X3t{^x)xU^y) + ^ i={*y)- (2-3) 

Here r counts the number of link operators contributing at each order. The coefficients c^'" ^ 
are the couplings between the operators S&, £. sitting at nearest-neighbor (NN) links l[ = (jc,-,y,-) 
in representation Mi. The effective action hence describes a network of link operators that are 
collected into (possibly disconnected) 'polymers' contributing with 'weight' c^j" ^ . One expects 
the 'weights' or couplings to decrease as the dimensions of the involved representations and inter- 
link distances increase. In a strong coupling (small j8) expansion truncated at ') one has 

r <k and the additional restriction \&i | H h \& r \ < k with |^| = pi for a given representation 

3% of SU (AO with Dynkin labels [p u ...,p N ]. 

Effective models for SU(3) Yang-Mills 

To lowest order & (fi Nt ) one finds the universal effective action 

Seff = Clo£S 10 ,te,) =hZ(^^; + ^*^y). (2-4) 
(xy) (xy) 

Our truncated model to order ^(/3 2A,t ) and with nearest neighbor interactions reads as 



Seff = Ai ZiXwWXoiW + cc.)+A 2 £(l6oWtoW + cc.) 

[xy) [xy) 

+ A3l(Zio(^)Zoi(^)+cc.) 2 . 

(xy) 

For a discussion of effective SU (3) Polyakov loop models see [||, §]. 

2 



(2.5) 



IMC and Demon Methods for Effective Polyakov Loop Models ofSU(N)-YM 



Christian Wozar 



3. Inverse Monte-Carlo - the basics 

The inverse Monte-Carlo (IMC) method [|j] allows to determine (effective) actions from given 
configurations. In our case, these are Polyakov loops obtained from gauge configurations generated 
with the Wilson action. Via IMC we determine the couplings of truncated effective actions which 
(ideally) would give rise to the same distribution of Polyakov loop configurations. 

The IMC procedure is based on an ansatz for the effective action of the type S e s = £;AjS,-. 
Translational invariance of the reduced Haar measure leads to Schwinger-Dyson (SD) equations 
[ph. They constitute an overdetermined linear system for the effective couplings A,- which may be 



solved by least-square methods (see m, K 10]). A second way to determine the couplings A; is the 



demon method which was successfully applied to SU (2) YM in [1 1 1. 



4. Microcanonical demon method 

Based on the large volume relation between microcanonical and canonical ensemble in statistical 
physics additional degrees of freedom ( "demons ") are used to simulate the effective theory together 



with the demons at a fixed total energy/action [12]. The demon is used as a "thermometer" to 



measure the coupling of the corresponding part of the effective action. 

To simulate a microcanonical system with action S[£P] = L; A/S,-^] we transform the canon- 
ical measure to the microcanonical one, 

p[&,E D ] oc exp(-£A,.(S ; -[^] +4)) — > 8[Si[&] + E\, - J&J. (4.1) 

i 

Each demons' energy is distributed according to p{E\ ) ) oc exp(— A^E^), A, depending on 
{Si}. The constraint E^ G [—E l ,E l ], E l Q < oo leads to an invertible relation (Ej^) = /i(A,). Couplings 
are obtained via the demon method in the following way: 

1. Simulate the microscopic (full YM) system without additional demons. 

2. Reduce the system for a chosen configuration to a Polyakov loop configuration. 

3. Perform a microcanonical simulation of the reduced Polyakov loop system with coupled 
demons. As discussed below the thermalization procedure should be handled cautious. 

4. The mean energy of the demons is directly related to the couplings of the effective theory. 



Tuning the microcanonical demon 

In the microcanonical method one YM configuration is reduced to a Polyakov loop configuration 
to start the microcanonical run. Therefore the method is highly sensitive to the chosen starting 
configuration. Thus, the effect of choosing specific (well thermalized) configurations with 5 ; | confi 
in the vicinity of (S;)ym * s analyzed below. 

For small E' thermalization problems of the demon arise due to the small acceptance in the 
update procedure. If E' is too large the demon is able to (and generically does) take away much en- 
ergy from the effective system. Thus, configurations within the microcanonical ensemble become 
independent of the starting configuration after the reduction step. These problems are circumvented 
in the following way: 

1. Choose a large energy range [—E' ,E' ] of the order ^(|(5') YM |). 

2. Reduce the YM configuration to the Polyakov loop configuration Co- 

3. For a few times (10 in our case) perform microcanonical simulations with Co as input con- 
figuration in every run. The demons' start energies are given by the expectation value (E^} 
in the preceding run. 

4. The final run lasts for the same Monte-Carlo time as the preceding runs and is used to mea- 
sure (ft). 
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Finally the contact between microscopic and microcanonical system is based only on one 
configuration. Further improvements should be possible by using a canonical demon method with 
improved "thermal contact" to the microscopic system. 

5. Canonical demon method 



In order to use the full statistics of the microscopic system we apply the following algorithm []13j] : 

1. Simulate the microscopic system according to e~ Sw until thermalization. 

2. Perform the reduction of the microscopic system to the effective system. 

3. Perform N m i cm microcanonical updates of the joined system of effective model and demon 
energies. These updates do not change the total energy S' +E^. 

4. Freeze the demon system and update the microscopic fields up to a new independent config- 
uration. After that proceed again with step ^[ 

To deal with thermalization effects of the demons' energies we begin the measurement after /^thermal 
microscopic configurations with a suitably chosen iVthermai- 

6. Observables for SU(3) 

We discuss the YM theory on a x_/V t -lattice. The Polyakov loop & x is measured in terms of its 
lattice average, 



1 



■ Nl 



(6.1) 



Since we deal with phases where the traced 
Polyakov loop is located halfway between the SU (3) l 
center elements we project the value of the traced j m p 
Polyakov loop onto the nearest Z^-axis and define a 
rotated Polyakov loop by (see Fig. [T]) 



ReP 
±ReP+^ImP 
JrReP-^ImP 



PE J?" 



-1 



-2 



(6.2) 




o 1 
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Figure 1: Fundamental domains of SU(3). 



7. The SU(3) YM phase diagram with adjoint potential 

An adjoint particle with mass M and spin s at temperature T leads to an effective potential Jl4| 

~(2s+l)M 2 T 2 



-K 2 {M/T)xn 



Thxn 



h<0. 



(7.1) 



Nevertheless topological excitations allow for a positive h and we therefore study a lattice action 

5 = P I ( l ~ ^ Re tr tfn) +#X>iiW ( 7 - 2 ) 

with standard SU (3) Wilson action and adjoint potential with unconstrained parameter H. 

Simulations of this system near the confinement-deconfinement phase transition on a 12 3 x2 
lattice with varying j3 and H show the phase diagram (Fig. |2|) in terms of the rotated Polyakov loop 
P ro t- At H = the well-known undirected and center-directed Polyakov loop structures related to 
(de)confinement appear. In the lower half-plane an additional structure arises where the Polyakov 
loop points into "anti-center" direction. 



The additional potential term is already contained in the effective model ( |2.5[ ). We therefore 
not only analyze the (de)confmement transition, but also look for a sensible analysis of the anti- 
center phase. 
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7.1 The confinement-deconfinement transition 

After simulating the (de)confinement phase 
transition (the upper black curve in Fig. HJ) 
we used IMC with the SD equations as 
well as the (micro)canonical demon me- 
thod for deriving couplings of the trun- 
cated effective theory. The programming 
codes were checked by simulating effec- 
tive theories with fixed couplings and re- 
producing them consistently with the SD (^rot 
equations and demon methods. 

The computed couplings correspond- 
ing to one YM coupling j8 are then used 
to simulate the associated effective the- 
ories with a Metropolis algorithm. The 
resulting expectation values of the rotated 
Polyakov loop are given in Fig. ||] (left 
panel). Here the SD equations fail to re- 
produce the phase transition point whereas the demon methods reproduce (P mt ) near the phase 
transition showing a better behavior than the SD method in the vicinity of the critical coupling. 




Figure 2: Phase diagram of SU (3) YM with adjoint potential 



according to Eq. (7.2) 
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Figure 3: Expectation values of the full theory compared to values produced with effective models after 
applying IMC methods. Left: (De)confinement transition at H = 0. Right: Anti-center phase at /3 = 6.05. 



7.2 The anti-center phase 

At p = 6.05 (vertical black curve in Fig. ^) we analyzed the anti-center phase and compared the 
resulting expectation values of P ro t- Whereas the methods show a quantitative difference in render- 
ing the critical point near the (de)confinement transition the expectation values of P rot indicate a 
qualitative difference in the present case (Fig. ||, right panel). While the SD equations produce a 
smooth behavior with high accuracy deep in the deconfined phase they fail to reproduce the anti- 
center phase completely. In contrast demon methods are sensible to the full phase structure. This 
behavior in the anti-center phase was analyzed with 750 microcanonical runs at j8 = 6.05, H = 0.2 
using different randomly chosen starts. The resulting couplings are plotted in the three coupling 
phase diagram of the effective theory which shows a symmetric, center directed and anti-center 
directed phase (see Fig. |j). Even far away from a phase transition in the microscopic theory the 
corresponding effective theory can be located in the vicinity of a phase transition of the effective 
model. 
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8. Using thermalized configurations for the microcanonical demon 

In the microcanonical demon method the starting configuration can be chosen randomly from the 
full YM ensemble. In contrast we can take configurations with 5,-| confi ~ (S;')ym ("well thermal- 
ized"). In the deconfined phase at j3 = 6.05, H = the projection of the three couplings to the 
A1-A2 plane comparing the different starts is shown in Fig. |5] (left panel). 

The extent of the coupling distribution 
shrinks for the well thermalized configurations 
compared to the randomly chosen configura- 
tions. Additionally the couplings of the well 
thermalized configurations correspond to al- 
most the same rotated Polyakov loop whereas 
the couplings derived from the randomly cho- 
sen configurations show a much broader dis- 
tribution of the corresponding rotated Polya- 
kov loops. 
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9. Thermalization 

effects of the canonical demon 



Figure 4: Phase boundaries of the three coupling effec- 
tive model and couplings derived from the microcanon- 
ical method at J3 = 6.05 and H = 0.2. 



One parameter of the canonical demon method is the number of microcanonical sweeps per micro- 
scopic configuration N m i cro . We analyzed the distribution of the couplings in order to read off any 
thermalization effects (Fig. |5[ right panel). The microcanonical distribution serves as a reference. 

Obviously the thermalization of the effective model has to be taken into account. With a 
small A'micro not enough time is spent to thermalize the effective system completely. The measured 
couplings refer to a non-equilibrium state of the effective model. Only in the large A^ m i cr o limit the 
couplings describe the thermal equilibrium of the effective model. This behavior is given by the 
fact that a configuration taken from a thermalized ensemble of the microscopic (YM) system is not 
necessarily a representative of an equilibrium state of the effective theory Jl5|]. Additionally the 
couplings obtained via Schwinger-Dyson equations do not correspond to the ones computed with 
the canonical demon method. 
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Figure 5: Left: Couplings obtained via the microcanonical demon method with randomly chosen and well 
thermalized configurations. Right: Thermalization effects due to different iV m i cro of the canonical demon 
method. 



10. Outlook to SU(4) YM 

For the finite temperature phase transition of the SU (4) YM theory on a 6 3 x2 lattice we applied 
the IMC method with the canonical demon method. The effective model is a generalization of the 
three coupling model for SU (3) YM, 
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Seff = M £ (ZiooWZooi i&y) +cc.) + A 2 £ (Zoio(^)Zoio(^y) + c.c.) 

W to) (1Q1) 

+ A 3 £(Z20o(^koo2(^) + c.c.) + A 4 £tooo(^kooi(^)+c.c.) 2 . 
to) to) 
Even for St/ (4) the demon method is a robust way of obtaining couplings of effective Polyakov 
loop models describing the phase transition (see Fig. 0). 
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11. Conclusions 

We studied and compared two ways of obtain- 
ing couplings for effective Polyakov loop mod- 
els which are related to SU (N) YM theories. 
The application of the inverse Monte-Carlo me- 
thod with SD equations to the St/ (3) YM case 
leads to stable results only far away from the 
phase transition. In the vicinity of the phase 
transition demon methods lead to a much bet- 
ter sampling of expectation values of the Polya- 
kov loop. We tried to reproduce the anti-center 
phase of a model with standard Wilson action 
and adjoint Polyakov loop potential by SD and 
demon methods. The SD method fails to reco- 
ver the phase structure while demon methods are favorable even in this case. The demon method 
can be generalized straightforwardly to the case of SU (4) YM leading to robust results in the 
vicinity of the finite temperature phase transition. Combining our experiences with both methods 
SD equations are less efficient than demon methods near first order phase transitions although SD 
equations have proven to be very useful for the second order transition in SU (2) YM 

When using demon methods much care has to be taken of different thermalization effects. 
Firstly with microcanonical demons the way of choosing microscopic configurations influences 
the derived couplings. Secondly when using the canonical demon method thermalization effects 
(which cannot be cured as discussed in Jl5|]) must be taken into account. 
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Figure 6: Finite temperature phase transition of 
St/ (4) YM on a 6 3 x2 lattice. 
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